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Abstract. We study the first passage dynamics of an ageing stochastic process in the 
continuous time random walk (GTRW) framework. In such GTRW processes the test 
particle performs a random walk, in which successive steps are separated by random 
waiting times distributed in terms of the waiting time probability density function 
~ (0 < a < 2). An ageing stochastic process is defined by the explicit 

dependence of its dynamic quantities on the ageing time ta, the time elapsed between 
its preparation and the start of the observation. Subdiffusive ageing CTRWs describe 
systems such as charge carriers in amorphous semiconductors, tracer dispersion in 
geological and biological systems, or the dynamics of blinking quantum dots. We derive 
the exact forms of the first passage time density for an ageing subdiffusive GTRW in 
the semi-infinite, confined, and biased case, finding different scaling regimes for weakly, 
intermediately, and strongly aged systems: these regimes, with different scaling laws, 
are also found when the scaling exponent is in the range 1 < a < 2, for sufficiently long 
ta- We compare our results with the ageing motion of a test particle in a quenched 
energy landscape. We test our theoretical results in the quenched landscape against 
simulations: only when the bias is strong enough, the correlations from returning to 
previously visited sites become insignificant and the results approach the aging GTRW 
results. With small or without bias, the ageing effects disappear and a change in 
the exponent compared to the case of a completely annealed landscape can be found, 
reflecting the build-up of correlations in the quenched landscape. 


1. Introduction 

Anomalous diffusion is characterised by the non-linear time dependence of the mean 
squared displacement (MSD) of the underlying stochastic process [H El [3], and it is 
widespread in the microscopic world. It can be detected in numerous systems, ranging 
form moving charge carriers in amorphous semiconductors BE], geological BE] and 
biological 0[9lEniEIlE2lE3lEl] systems, and similar features also occur in the dynamics 
of blinking quantum dots [15] and in laser cooling [T6] . 

There are several models describing subdiffusion [I1131E71EH1E9] and even ultraslow 
logarithmic diffusion [201 EB E2l EB El]- The continuous time random walk (CTRW) 
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ageing period ta measurement t 


Figure 1. Sketch of an ageing subdiffusive process. The observation is started after an 
initial ageing period of length ta, meaning that the origin of the observation time t lies 
in the point ta- When the scaling exponent is in the range 0 < a < 1 statistically longer 
and longer waiting times appear in the evolution of the process. In this non-stationary 
scenario physical observables explicitly depend on ta- 


is one of the best studied anomalous diffusion processes BE]. It exhibits subdiffusive 
behaviour 

(x^(t)) ~ (1) 

with 0 < a < 1 for power law waiting time distributions 

~ ( 2 ) 


of the waiting times r elapsing between any two jumps. Successive waiting times 
are statistically independent of each other. For 0 < a < 1 the distribution 
is characterised by an inhnite mean (r). Below we will also consider the case when 
1 < a < 2 and thus (r) < oo, yet the fluctuations around (r) diverge. Theoretical 
dynamic models that can be effectively analysed within the CTRW theory include 
dynamic maps [251 [26] and particles moving in random (or quenched) energy landscapes 
[311271 EH]. CTRWs exhibit several interesting features, such as (i) weak ergodicity 
breaking, the violation of the Boltzmann-Khinchin equality of the ensemble and lag time 
averages of physical observables [HI EHl [301 EH EH [33l EH [35], (h) population splitting, 
the existence of a growing immobile fraction of particles versus mobile particles [361EZ] , 
and (hi) ageing, the dependence of the time evolution of statistical properties on the time 
lag ta between preparation and start of the observation of the system [261 EH EH EH] , see 
hgure [TJ Even more signihcantly, the scaling of physical observables with time changes 
with the age, as can be seen, for instance, in the MSD, which turns from subdiffusive 
growth for weak ageing ta t to seemingly normal (linear) growth for strong ageing 

4 > t [HHEH], 


{X^{ta,t)) ~ 


Kat"", ta^t 
ta:>t 


(3) 


Some ageing effects persist also for waiting times distributed according to equation ([2]) 
with 1 < a < 2 [39l [iQ]. Evidence of ageing was found experimentally in biological 
systems mm as well as in amorphous semiconductors [12|, glasses [30], and in the 
blinking dynamics of quantum dots [15] • Note that ageing and ergodicity breaking 
also occurs in other stochastic processes such as Sinai diffusion 1221 , scaled Brownian 
motion [IH [H], processes with spatially heterogenous diffusivities [IH [16], as well as 
their combinations (generalised diffusion processes [U]) and many-particle systems [23] . 
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The first passage time, i. e., the time a process needs to reach or cross a certain 
threshold, is an important tool in the statistical analysis of stochastic processes. The 
hrst passage describes, for example, the arrival of a particle at an absorbing barrier or 
to the binding of a protein to its binding site. The scope of this paper is to describe the 
effect of aging on the hrst passage time density (FPTD) for a one dimensional CTRW 
process in a semi-inhnite domain (one absorbing boundary) and in a hnite domain (the 
particle starts between two absorbing boundaries), thereby extending and completing 
our previous work [IH]. In section [2l the FPTD for the unbiased cases is derived. The 
cases of 0 < a < 1 and 1 < a < 2 are considered separately. In section [3] we continue 
our work on the FPTD for biased diffusion j3H] and section 0] is devoted to the ageing of 
the FPTD in annealed and quenched random energy landscapes. The paper ends with 
a discussion of the relevance of ageing in hrst passage processes. 


2. First passage time density for ageing CTRW in semi-infinite and finite 
domains 


The general idea for the derivation of the hrst passage time density (FPTD) |39] is to 
use the survival probability of the random walker [50] in analogy to the non-aging case 
[5T| . The survival probability ^{t) is given by the integral of the probability density 
function Po^xfiafi) to hnd the particle at position x at the time t after the ageing 
period ta over the spatial domain D, in which the walker is conhned to move. 


^(s, u) 


Pd{x, s, u)dx. 


D 


(4) 


For simplicity in the following we use the double Laplace transform {ta, t ^ s,u) of the 
quantities of interest, i.e. 


Pd (x, s, u) 


— Ut—Stc 


Pd {x,ta,t) dtdta- 


0 0 


(5) 


The FPTD p(ta, t) is then given by the negative time derivative of S^{ta, t), or in double 
Laplace space as 

p{s,u) = - uy'{s,u). (6) 

s 

The propagator for the unconhned aging CTRW with initial condition P {x, ta,0) = 6{x) 
is well known. In double Laplace representation it reads [3H] 


P{x, s, u) = <h(s, u)5{x) + h{s, u)Pmw{x, u). 


(7) 


The time t we have to wait for the hrst jump to occur is statistically drawn from the 
distribution 


1 fi’is) — fi{u) 


l-fi{s) 


h{s, u) 


u — s 


( 8 ) 
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of the forward waiting (or recurrence) time, whose statistic is different from that of 4>{t) 
[39] . In equation ([7]) the PDF h {ta, t) is convoluted with the Montroll-Weiss propagator 
PMw{x,t) [52] for the non-ageing CTRW. The first term on the right hand side of 
equation ([7]) is the probability of not having made a jump between ta and t, 


N 1 h{s,u) 

su u 


(9) 


The probability density functions in the semi-inhnite domain (with and without force) 
and in the finite domain all have the following general form 


Pd{x, s, u) = h{s, u)Pd{x, u) -|- <I)(s, u)6{x). (10) 

Here P£,{x,u) is the known non-ageing solution in the respective cases [SniEI]. The 
combination of equations (Hj), (Ej), (I9|) and (ITOl) yields a general expression for the ageing 
FPTD, 

p(s, u) = -m<F(s, u) — uh(s, u)y'iu) 

s 

= h{s,u)\i — uS^iu)] = h{s,u)p{u) ( 11 ) 

This result implies, that the aging FPTD is the convolution (the Laplace space product) 
of the forward waiting time density h {ta, t) and the non-aging FPTD. 


2.1. Semi-infinite domain 

The FPTD of a non-aged CTRW on the semi-axis, starting at x = 0 and with the 
absorbing boundary located at x = xq, is proportional to a one-sided Levy stable law 
(Levy-Smirnov law) [53] 




( 12 ) 


which is defined via its Laplace transform in the form 


p(u) = exp 




(13) 


The forward waiting time distribution has different expressions for the two ranges of a 
we consider, respectively. For 0 < a < 1 one hnds the form [371 139] 


h{ta,t) 


sin (cTTr) f ta\°' 1 


(14) 


which is valid for all ageing times. For 1 < a < 2 the mean waiting time is finite 
((r) < oo) but the fluctuations are unbound ((r^) — )■ cxo) [5l]. In this case the 
asymptotic forms are known, which are h{t) = 1 for ^ 0, while for ta ^ oo [371 EH] 
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This is not sufficient to find all regimes of the FPTD. We use the full form ([H]) in the 
Laplace domain instead, in combination with the Laplace space waiting time density 


fi{u) ~ 1 + (f) M + r“r[l — «]«“, 


(16) 


valid in the limit m —?■ 0. For a > 2, ageing effects are formally present, but decay 
too fast to be detectable [37], compare similar transient ageing and weak ergodicity 
breaking effects in asymptotically ergodic systems [SSlEnilST]. We first deal with the 
case 0 < a < 1 and come back to the second case below. The central result for the 
aging FPTD is the Laplace convolution in t hhiQ 


P{ta,t) = 


sm air 


71 


ta t 




Xn 


£a 


K„ 


(17) 


which follows directly from equations flTT]) . flT^ . flT4|) and the convolution theorem of 
the Laplace transformation [5H]. For practical reasons one can derive a result in form 
of a series expansion of equation flTTl) in powers of f, using the single Laplace transform 
of h(ta,t) with respect to t, 


h{ta,u) = exp{tau)T (afiau)'^ (a 


(18) 


A closed form solution of the inverse of the full Laplace space result cannot be found. 
The incomplete F-function can be expanded in a series 


F {a,tau) = F(q;) 


exp {-tau) ^ 


{taU, 


.n+a 


—^ F(1 + a + n) 

n=0 ^ ' 


(19) 


Similarly, the series for the non-aging FPTD is 

Xq 




m=0 




u 


maf2 


( 20 ) 


With these series, equation flTT]) turns into 


p{ta, u) = exp (taU) V (-^ 

m=0 V V « 


U 


mal2 


- E 


-j-n+a 


^ m!P(l + n + a) 

771,71 = 0 ^ ' 


3^0 


U 


n+Q;+mQ;/2 


( 21 ) 


The inverse Laplace transform of this expression can be obtained term-by-term. The 
exponential function yields a shifted Laplace transform u ^ t + ta oi the remainder and 
the powers of u are inverted using the Tauberian theorem [59] 




t 


- 1-6 


r(-6) 


( 22 ) 


I I.e., f{t) O g{t) is defined as / g(T)f{t — T)dr. 

0 
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This formula for the inverse Laplace transform ^ ^ is valid for u —)■ 0 and t ^ oo, and 
thus we find the aging FPTD 


Pi^a, t) — ^{ta + t) + 

m=l 


f _xq + 

V y/T^J m\T{—ma/2) 






m\T(l + n + a) 

n,m=0 ^ ^ 


Xq 


^ —1—n—Q;(H-m/2) 


y/K^J r [-n - «(! + m/2)] ‘ 


(23) 


The asymptotic behaviour for f —)■ cxo, meaning the dominating powers of t, are found 
from comparison of the terms with the smallest exponents in t, assuming that we are in 
the weak ageing regime tta- We find two different regimes with a crossover at 


t* = t 


2 

a 


^y/KaT (1 — a/2) sin (vra) “ 

XQTra 


(24) 


At t*, the ageing effects become negligible and the final non-ageing regime is reached. 
Note that the crossover time t* depends on the initial condition xq. In the strong ageing 
limit ta t one has to employ the asymptotic expansion of the incomplete T-function 
for uta —)■ oo, which is found using the integral dehnition and partial integration — this 
result cannot be found from the Taylor series. The dominating term is: 


T (a, tau) ^ {utaY ^ exp {-uta), 


(25) 


for uta —)■ cxD. In this limit the forward waiting time is dominating the rest of the hrst 
passage process, and consequently, the inverse Laplace transform of this asymptotic 
FPTD yields to lowest order 


P{ia,t) 


sm TTCK 


TT 




(26) 


This is the just the limit of h{ta,t) for ta 3> t. We have thus identihed three regimes, 
separated by the crossover times ta and t*, namely 




P{ta,t) ~ < 


sin (7rQ()^„_i^_„ 

TT 

TT 

_ ^0^ _^-l-a/2 

2r{l-a/2)y/KT 


ta'>t 

ta^t<^t* 

t* 


(27) 


The hrst exponent of t is —a, the steepest exponent —1 — a can be found in the 
intermediate time regime. In the asymptotic regime one hnds the intermediate exponent 
— 1 — q:/2 equivalent to the non-ageing case [2l|50]. Since the crossover time t* decreases 
with the initial distance Xq of the diffusing particle from the absorbing boundary, we 
hnd that the hnal regime takes over faster for larger xq . As a result, the probability 
density decays slower than for short distances, which is intuitively expected. The three 
regimes fl271) can be clearly distinguished in simulations, as shown in hgure [2l One can 
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Figure 2. FPTD in the semi-infinite case (blue line). One can see the three scaling 
regimes of equation (l?7)) . which are separated by ta = 10® and t* Ri 3.15x 10^° indicated 
by the vertical lines. The parameters are a = 5/7, xq = 1 and Ka = 1/5, green: ~ 

(1 "C ta), red: ~ (ta "C t t*), orange: ~ f-i-a /2 ^ dotted line is 

the non-ageing FPTD equation (fT^ . 


see that the non-ageing power law is only reached after extremely long times. On short 
or intermediate time scales, the ageing effects are dominant and if the corresponding 
parameters allow both power laws persist long enough to be in principle detectable in 
experiments. 

In the case of 1 < a < 2, for which the mean waiting time (r) is hnite, one hnds the 
long time limit of the FPTD in the non-ageing case to be the same Sparre-Andersen 
law as for Brownian diffusion (hgure |3]) , 

_^ 0 _^- 3/2 

2\/ Kir 




For the ageing case, one uses the general form (IH)) of the forward waiting time and the 
waiting time density (1T6|) . The ageing FPTD (cf. equation (ITT]) ) is analysed in a similar 
way as for 0 < a < 1. As a hrst result, three power laws can be distinguished for 
1 < cc < 3/2, which are separated by the crossover times ta and 


* 




2/(2a-l) 


t 


2ax (t) 


(29) 
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which is again the time that separates the ageing regimes from the non-ageing one. The 
power law behaviour in those regimes is, respectively, 


P{ta,t) ~ < 





{t) a 


j^—l — OL 


_^0_^_3/2 

2-\/ i^TT 


tai^t 

ta<^t<^t* 

t* 


(30) 


The first regime is the one, in which the forward waiting time density dominates the 
FPTD (cf. equation flT^ i. The last power law is that of the Brownian limit. In the time 
range f the system has evolved for so long that the fluctuations of the aged waiting 
times around (r) become negligible. Only in the crossover regime the dominating term 
depends on the ageing time explicitly. The three regimes are shown in figure El where the 
middle regime is not indicated by a power law line as the picture would be overloaded. 
For a > 3/2, these regimes collapse to a single non-ageing regime (cf. figure E]), which 
is just the Brownian one (equation (l28|) b We will see in sectionEl that the introduction 
of a constant bias will change the long time behaviour severely. 


2.2. Bounded domain 


In a bounded domain [—a, a], we find the solution of the non-ageing FPTD using a 
trigonometric basis for expansion that vanishes at the boundaries x = ±a. This is 
equivalent to an infinite application of the method of images [59]. For the symmetric 
initial condition P{x,ta,t = 0) = S{x), the solution reads 


Pd{x, s,u) 


1 

a 


Re 

m=0 



f .{2m -\- l)7r 

I ^ 

V 2a 



P k 


{2m -I- 1)71 
2a 



(31) 


where Re {■} is the real part of the argument. To solve the asymmetric initial condition 
P{x,ta,t = 0) = 6{x — Xq), —a < Xq < a we shift the sum to the coordinate system 
x' = X + xo and introduce the factors am, so that 


1 oo . 

Pd{x', s,u) = — Re < am exp 

^ m=0 ^ 


.(2m -I- Ijvr 
* ^ 


[x 


Xo) 


p k = 


{2m + l)7r 


(32) 

This solution is in fact the propagator for the symmetric initial condition P{x,ta,t = 
0) = I [5(a; — Xq) + 6{x + Xq)], but the resulting FPTD is the same as for the asymmetric 
case. A way to find the propagator for the asymmetric initial condition is mentioned in 
Ref. [19]. With equation (E2]), we have to match the 6 initial condition (at x' = 0, due 
to the coordinate shift) and thus we identify 


.(2m -I- l)7r 



am = exp 


2a 


(33) 
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Figure 3. Numerical Laplace inversion of equation ([5]) multiplied with the Brownian 
FPTD p{u) = exp for a = 1.1 (blue, dashed line, top) and a = 1.9 

(red, dashed-dotted line, bottom) in the strong ageing limit (the Laplace variable is 
chosen as s = 10“®, which roughly corresponds to an ageing time ta ~ 10®). The first 
case with a = 1.1 is compared to the forward waiting time (1151) (green), the second 
case with a = 1.9 to the Brownian FPTD (orange), see the top right corner: same 
plots, larger scale: all three regimes (cf. equation (1301) 1 are shown, the kink is due to 
numerical errors in the Laplace inversion. The function plotted is p(t, s) with fixed 
s, which behaves similar to p{t,ta), since in the hrst and in the final time regime the 
value does not depend on the second variable. We only have one less numerical inverse 
Laplace transformation and thus one less source of numerical errors. 


The result in the original coordinates is now 


Pd{x,s,u) = 


OO ^ 

m=0 ^ 


.(2m + l)7r 

exp ( ^- — - {x + xo) 


(2m+l)x 
2a ’ 


(34) 


We replace the ageing unconhned propagator using equation ([7]) 
h{s, u] 


Pd{x,s,u) = 


OO ^ 

J^Rej 

m=0 ^ 


.(2m + l)7r 

exp ( t - — -(x — xq) 


<^{s,u] 


OO ^ 

E R* { 

m=0 ^ 


exp i 


2a 

.(2m + l)7r 


MW 


k = 


(2m + l)7r 
2a 


,u 


2a 


{x - Xo) 


(35) 


In the hrst line, we just have the product h(s,u)P£)(x,u) on the right hand side. The 
sum of the exponentials can be transformed into a sum of ^-functions, using the Poisson 
summation formula, and 


OO —1 

exp (iTTnx) = exp(i7ma:). 

n=0 


n=—oo 


( 36 ) 
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f .{2m+ 1)71 , ^ 

exp 7 ---(x - xo) = 


m=0 




2a 


} = ^ £ Re|exp('z 

m —_ r>ir> ^ '' 


(2m + l)7r 


2a 


{X - Xq) 


m=—OQ 

OO 


= a 


S{x — Xq — 2a). 


(37) 


m=—OO 


Only the peak for m = 0 is within the bonndaries of the domain and thus Pd{x,s,u) 
has the form of equation (fTOj) . with a 6 initial condition at x = Xq instead of x = 0. But 
since the integral for the survival probability is taken over the whole domain, this does 
in fact not make a difference for the FPTD. In Laplace space, we obtain jSO] 

p(ta, u) = r (a, tau) exp 

4(-l)"*exp {i 


X 


l-J^Re 


u 


m=0 


(2m + l)7r 


(2m+l)27r2 




(38) 


To hnd the solution of the non-ageing part in the form of a power series, we expand the 
last fraction in equation fIMl) about u = 0, 


u 


Aa u 


2„,a 


Aa u' 


2.,a 




I ( 2 m+l )2 7r2 iL (2m -|- li^TT^ ^ 

' 4^2 / r)=i 



_JL„(2m+ 1)%2 

n+1 


n=l 


Changing the order of summation, we get 

p[ta, u) ^ r (a, taU) exp (^„^^)^"^(tt) 


(39) 


(40) 


X 


1 + E 


n=l 


2„,a 


a U 




J^^Re 

m=0 


22’^+2(-l)™exp 


(2m l)2n+l7]-2n+l 



As the series in m equals unity for n = 0, 

4(-l)-exp(z(^^Xo)l 4 


^Re 

m=0 


(2m -|- l)7r 


= —Re <1 arctan ( 


71 

= 1 , 


2a 


(41) 


for —a < Xq < a, we can write the 1 in equation (1401) as the zeroth element of the series 
in n. Furthermore, the series in m can be written as a Lerch function ^ |60] : 


^Re 

m=0 


22’^+2(-l)"^exp 


(2m -|- l)2n-|-1^2n+l 


I 


y^2nH-l I 


Xq 


Re < exp ( iTT— )R 


2a 


0 - + X 0 1 

-, 2n 4- 1, - 

2a ’ ’2 


(42) 
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so that the FPTD as a series in powers of u°‘ reads 
p(ta,u) = r (a,4M)exp 


11 


E' 


X > u" 

n=0 


a 


TT 


2nH-l 


Re |exp 


a + xo ^ , , 1 

-, 2n + 1, - 

2a ’ ’2 


(43) 


Using also equation (IT^ . the Laplace inversion u ^ t yields 

j.n+a / 2 \ ™ 

p(t.. () = P(K + () - E r(i;„+a) {-irj 

r) ryi / \ / 


X 




I 


Re < exp 


r [—n — a(l + m)] ( 


tTVXo 

2a 




a + xo 1 

-, 2?7, + 1, - 

2a 2 


(44) 


In the symmetric case (xq = 0), one can simplify this due to the relation between the 
Lerch function and the even Euler numbers [60] 


y^27TT.“|~l 


i? (1/2, 2m + 1,1/2) = 


IE. 


2m 


(2m) 


(45) 


It should be noted, that the symmetric result is more intuitively found when the 
coefficients are written as an inhnite series of cosines instead of Lerch functions: 


7j-2n+l 


Re < exp 


tTlXo 

2a 




a + xo 1 

-, 277 , + 1, - 

2a 2 




E 


cos (vm + ^ [1 + 2n]) 


l\2n+l 


(46) 


m=0 (™ + I)' 

Setting xq to zero, we are left with the cosines of mr, which is just (—I)"". This alternating 
series is exactly the series for the Euler numbers. 


EonI — 


22’"+2(2n)! 


E 


-1)’ 


^2n+l (-2^7^ + l)2n+r 

m=0 ^ ^ 


(47) 


The behaviour of the FPTD for short ageing times is 
p{ta,t) ~ 

F(—ck)F(1 + a) 

''inxo 


E„7r3p(-a)' 


Re < exp 


2a 




a + xo 
2a 


,3,1/2 


or in the symmetric case 


P(—a)P(l + Q:) 2 E'q,F(— a) 


(48) 


(49) 


If ta is sufficiently large, the second term in the square brackets, and thus also any 
asymmetry effect, becomes irrelevant. In the strong ageing limit ta ^ t (cf. figure Hj) 
the behaviour of the FPTD is dominated again by the forward waiting time, and thus 
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'H' 


'H' 








t 


Figure 4. Simulation (blue squares) and exact FPTD (red line) of an ageing CTRW 
in a box (green (lower) line: result without ageing); the power laws (~ t~°‘ for t ta 
and ~ for t ^ ta) are indicated by black lines, the crossover at ta is indicated 

by a vertical black line (ta = 2.5 x 10^, a = 4/7, = 1/2, a = 1, and ~ 10® runs). 


the result is the same as in the semi-inhnite case (cf. equation fl26|) ). The power laws 
in the long time limit are thus (cf. Ref. |1H]) 



(50) 


One can see that a nonzero ageing time will yield an offset in the limit ta t, and 
different from the semi-inhnite case, the non-ageing and ageing curves do not collapse 
onto the same curve (cf. hgures|l]and|2]). In the case 1 < ct < 2, the non-ageing FPTD is 
approximately an exponential distribution, which is found in the case of normal diffusion 
(similar to the results derived in [6T]), and yields a hnite mean hrst passage time. This 
means, that, after the hrst step is made, the particle will be absorbed very fast. The 
long time behaviour is thus dominated by the forward waiting time density for long 
ageing times (equation flTbll ) and an exponential decay in the case of a weakly aged 
CTRW. 

3. First passage density for biased ageing CTRW 

We now introduce a constant bias F acting on the random walker. This problem is 
handled using the method of images with corrections due to the external held, or by 
direct solution of the advection-dihusion equation ^91 |5T] . The propagator on the semi 
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axis reads 


Pd{x, s, u) = Pf{,x, s, u) — exp 


Fxr 


Pf{x — 2xo, s, u). 


(51) 


In this, the value of the Boltzmann constant is set to unity. The function Pf{x, s,u) is 
the unconfined ageing propagator under a constant bias, 


Pf{x, s, u) = <h(s, u)S{x) + h{s, u)Pf{x, u), 
and Pf{x,u) is its non-ageing equivalent 


Pf(x,u) = - , exp 

T^1 + A{ut(A)» 


F (^x — \/l + 4(Mrb))“|x| j 


2T 


(52) 


(53) 


with the scaling time /[F‘^Ka). From the additive form of equation (15T1) one 

sees immediately, that equation flT0|) remains valid. Using [51] 


p{u) = exp 





(54) 


one consequently finds the ageing FPTD via equation flTT]) . We can write the non-ageing 
result as an infinite convolution of one-sided Levy stable laws or as a series in powers of 
t. In the first case, we expand the square root in the exponent using [62] 


(l + (r)“-i 


•U r(a) . 
^na-m ■ 


x\ < 1. 


(55) 


An exponential to the power of a sum of terms can be written as a product of 
exponentials, and we thus find 


pip) = JJ exp 

k=l 



r(3/2) 
r( 3/2 - k)k\ 



(56) 


The exponential function is the Laplace pair of the one-sided Levy stable function 
according to equation ffT2]) and thus, the non-ageing FPTD is an infinite convolution of 
Levy stable functions 


pit) = (^Pi(y,k)ikaiPi<y,k)t) (57) 

k=l 


with 


T>ia, k) 


P(3/2) 


(58) 


§ We define f{t) as f{t) 0 f(t) 0 f(t )... 0 /(t). 
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The ageing FPTD is just the convolution with the forward waiting time. For the power 
series, we hrst expand the non-ageing FPTD 


p{u) = exp 


/ Fxc 


V 2T 


72=0 


Fx( 


2T / n\ 


rtrt\ L J 


(59) 


Then we use the series expansion (15^ again, only this time with the exponent depending 
on n. Change of order of summations yields 


p{u) = exp 


Fxq 

2T 


k 1 




k=0 


72=0 


Fxq 
' 2T 


1 r(^) 

n\ P - k) 


(60) 


The sum in n can be turned into a combination of hypergeometric functions. First 
we separate the even and odd numbered terms into two sums and change the index of 
summation. 


E 


(#)” r(!=±^) 


n=0 


n\ 


E 




2m 


F (m -|- 1) 


-E 


(^) 


2772 + 1 


F (m + I) 


i! r - k) ^^r(2m +l)r(m +1 - k) r(2m + 2) F (m + § - k)' 

(61) 

Then, using the duplication formula 


r(.- + i)_ r(i) 


r(2z) 

we expand the fractions 

(^r 


22-‘r(s)’ 


(62) 


n! 


E 

72=0 

oo 

E 


772=0 

_ 

AT 


Fxr 


2m -I 


iT 


m!r(i -k) r (1) 


T {m + 1 — k) m!F(l — fc) F (m -|- |) 


E 


777=0 


Fxr 


2m 


r(l 


0 r (y 


4T y (m -I- I — /c) F (I — /c) ml 
Finally, introducing Pochhammer symbols (x)(„) = P(x -|- n)/F{x), we hnd 

(^)“ r(!!±2) 


(63) 


E 


71=0 


n\ r (5±2 - k) 


r(i/ 2 ) 


rTiE 


OC) / j \ 2772 

Fxn\ 1 


P(l/2)P(l-fc)f-^V4T 




Fxq P(1/2) 

4T P(3/2-fc) \ AT 


E 


f FxqV"^ 1 




ml (1 - /c)(^) (l/2)(^) 
1 


ml (3/2 — k)^ 


= P 


1F2 ( 1 ; -, 1 — fc ; 


Fxq 

AT 


Fxq 

AT 


3^1 


Fxq 

AT 


(64) 
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Figure 5. Simulation of the FPTD of a biased ageing CTRW for three different 
anomalous exponents a = 0.4 (red, circles), a = 0.6 (blue, squares) and a = 0.9 
(black, stars). The parameters are F = 1, T = 1, ta = 2500, xo = 5 and Ka = 1/2. 
The lines are the theoretical results (p51) . the short black slopes indicate the power 
laws (1551) . 


in which the pFg are regularised hypergeometric functions |60] . The ageing FPTD reads 


pit, to) = p{t + to) + exp 


Fxq 

2T 


E 


^n+CK 


^-l-n-a{bm+l) / \ ^ 

r(l + n + «) T (—n — a\bm + 1]) \F‘^Ka) ml 


X 


1^2 I 1;-,1 - m ; 


Fxq 

AT 


Fxq 

AT 




3 

2 “ 


Fxq 

AT 


with the long time limits (cf. |1H]) 


(65) 


sm(7ra) 1 „ 


p{ta,t) = 


TT 




aTxn 


FK^T{1 - a) 


sm vra 


TT 


ta > t 


t ^ ta^t 


( 66 ) 


We again find an offset for ta t, which separates the ageing and the non-ageing curve. 
In figure [5] we simulated a biased hrst passage process for different anomalous exponents. 
We hnd that the power law behaviour in the strong and weak ageing regions agree well 
with our long time approximations fl66l) . For 1 < a < 2, we follow the arguments in 
section |2.2| and reference the FPTD is dominated by the forward waiting time 

density (IT^ for strong ageing and decays with a finite mean, and almost exponentially, 
for weak ageing. 
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4. Ageing first passage density in quenched and annealed energy landscapes 


We now consider random walks in a random energy landscape [271 EH EDI Ea] on a 
lattice. Each time the walker is updated it jumps to the left or to the right with equal 
probability in the unbiased case. At the new site, it is trapped with random energy 
E, which is drawn from an exponential distribution p{E) = exp {—E/Tg) [30] . 
The thermal (Kramers) escape from the energy traps dominates the particle motion; it 
yields a relation between the waiting time r and the trap energy E, which is given by 
the Arrhenius law [3] 



(67) 


where C is the inverse frequency of escape attempts. The distribution for r is thus a 
power law 

fiij) = ( 68 ) 


which yields subdiffusive behaviour, if a 
is 


= T/Tg < 1. The anomalous diffusion constant 


2r(l-a)C“’ 


(69) 


introducing here the lattice spacing g. If a constant bias is to be introduced, we just 
change the probability to jump to the left from one half to 1271 


1 

l+exp(Fg/r)' 


(70) 


We can have two different classes of energy landscapes. In the annealed random 
landscape, in every step a new energy is drawn from the exponential distribution. On 
the return to a previously visited site, the walker does not hnd the same energy as 
before. This case models exactly the subdiffusive CTRW and one hnds the subordination 
relation [21 [5T1 |6ll [65l [66] 

n ~ (71) 

between the number of jumps and the process time. This is the CTRW process 
considered in sections El and |3l In the case of a quenched energy landscape, where the 
energy of each site is hxed (correlations are introduced), the subordination is weakly 
broken [30] and one hnds [6^ 

(' 72 ) 

The change in behaviour is due to the multiple appearance of the same waiting times 
on revisiting specihc locations. A realisation of a quenched energy landscape is depicted 
in hgure El If the correlations in hnding the same waiting times repeatedly at the same 
position are sufficiently reduced, for instance by a large bias or in higher dimensions, one 
hnds the results converging to the annealed, renewal CTRW. We illustrate the change in 
the behaviour for CTRW and motion on a quenched landscape for diherent bias strength 
E in hgure [71 on the left one can see two diherent curves for small E, while coinciding 
results on the right appear for large E. 
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Figure 6. Realisation of an energy landscape with Tq = 1, T = 2/3 (top). Bottom: 
corresponding waiting times t(x) = exp {E{x)/T). 




t 
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Figure 7. FPTD for the quenched energy model in the semi-infinite domain: results 
of biased random walk simulations in an annealed (blue) and a quenched landscape 
(red) with different results for small bias F = 10“^ (left) and converging results for 
large bias F = 100 (right). The parameters are a = 1/2, ta = 100, T = 1/2, g = 10“^, 
N = 200, and C = 10“^. The data were averaged over ~ 2 x 10^ (left) and simlO® 
(right) runs. The straight blue line represents the theoretical result for the CTRW. 


For an explanation of the long time behavionr of the FPTD in a quenched energy 
landscape we use a result from the theory of Brownian random walks, according to 
which the number of jumps n is proportional to the time t. The asymptotic behaviour 
of the Brownian survival probability in the semi-inhnite domain is S^{t) ~ and 

thus we obtain 


S^{n) ~ n 


1/2 


(73) 
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Figure 8. FPTD in the semi infinite domain without external bias: results of annealed 
landscape simulations (left) and quenched landscape simulations (right) with (red dots) 
and without (green squares) ageing. The parameters are a = 5/7, fa = 2 x 10^, 
T = 5/7, g = 5 X 10“^, iV = 50 and C = 10“^, and averaging is taken over ^ 2 x 10"^ 
runs. The straight lines are the theoretical results ((23ll . the short black lines indicate 
the respective power laws and dill). The intermediate regime in the annealed 
landscape is too small to see the corresponding power law. 


Using equation (1721) . we expect for the quenched energy landscape that 


pit) 


dij 


(74) 


One can clearly see the predicted power law in our simulations in hgure [8] (right). 
Interestingly, the unbiased random walk in a quenched energy landscape exhibits no 
detectable ageing, unlike the random walk in a completely random energy landscape 
(hgure [HI graph on the left). This is due to the fact that the quenched energy landscape 
features only a comparatively small sample of random variables relative to the annealed 
landscape or a continuous time random walk drawn from the power law waiting time 
distribution ([2]). Ageing enters the system because of immobilisation of the walker for 
extremely long waiting times. Such extreme events are found in the large sample of 
random waiting times in a CTRW or in an annealed random landscape, but much fewer 
realisations occur in the hnite QEL. These ageing QEL trajectories are statistically 
negligible, as rehected in the FPTD observed here. 


5. Discussion 

We studied the hrst passage time distribution of an ageing CTRW in three important 
scenarios, namely the semi-inhnite case with and without external bias, and the conhned 
case. We presented the relevant details of the derivation of the exact form of the hrst 
passage time densities and extended our long time approximations from Ref. |18] to 
the asymmetric conhned case and to the range 1 < a < 2 of the waiting time scaling 
exponent. In the second part of the paper we show that a diherent scaling regime 
appears in a quenched energy landscape, which converges to the CTRW regime with an 
increasing bias towards the absorbing boundary. The efficiency of the hrst passage is 
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generally reduced by ageing. This effect is not only found in subdiffusive processes but 
is shown to exist even for waiting times of the form ([2]) with 1 < a < 3/2. Even though 
the corresponding CTRW shows normal behaviour, we can End an aged scaling regime 
in the strongly aged FPTD, which persists for many decades. For 3/2 < a < 2, we 
End no ageing in the unbiased semi-inEnite domain, but ageing exists in the bounded 
domain as well as for biased difiusion. The ageing efiects are so strong, that even a 
short ageing time is enough to overcome the efiects of asymmetry in a confined domain. 

The starting position of a particle in a box becomes (almost) irrelevant in the 
observable regimes. We find, that the process is accelerated in a quenched energy 
landscape, which counters the slowing down of ageing. In fact, ageing cannot be observed 
in a quenched energy landscapes without bias due to the finite sample of energy wells. 
The slope in the long time limit changes to —1 — a/(a + 1), but if a bias is introduced, 
the correlations are overcome and the CTRW regime is recovered excellently. 

First passage times are an important quantity for the evaluation of experiments 
and simulations, as they can be easily obtained and can also be used as diagnostic 
tools of stochastic processes [67]. Ageing is especially important in biophysics, 
condensed matter physics and geophysics exhibiting CTRW behaviour, where it 
corresponds to an immobilisation of particles in crowded cells, changing characteristics 
of electrical currents in amorphous semiconductors or decreased chemical release rates in 
groundwater aquifers. We extended the existing results on ageing FPTDs and examined 
the behaviour in quenched energy landscapes, which are a realistic model for many ID 
problems. Future work could deal with the problem of first passage in higher dimensions 
(crossing a plane or barrier). Also, it would be interesting to investigate the transition 
from biased non-ageing quenched energy landscapes to ageing CTRW motion with 
increasing bias. Another interesting related topic, the ageing Scher-Montroll transport, 
is and was examined theoretically [38] and also in time-of-fiight experiments [12]. The 
experiments in amorphous semiconductors promise excellent testing conditions for the 
CTRW model for charge carrier transport, as the age of the material can be reliably 
controlled. In groundwater systems mm the first passage times can be tracked over 
several decades, so that crossovers between different scaling regimes could be resolved 
with good resolution. In different systems, our results for the crossover times and 
scaling laws of different power law regimes could be used to determine the age ta or the 
stochastic nature of the underlying process. 
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